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Abstract 

We consider the natural super symmetry scenario in the framework of the iJ-parity conserving 
minimal supersymmetric standard model (called natural MSSM) and examine the observability 
of stop pair production at the LHC. We first scan the parameters of this scenario under various 
experimental constraints, including the SM-like Higgs boson mass, the indirect limits from precision 
electroweak data and B-decays. Then in the allowed parameter space we study the stop pair 
production at the LHC followed by the stop decay into a top quark plus a lightest neutralino or 
into a bottom quark plus a chargino. From detailed Monte Carlo simulations of the signals and 
backgrounds, we find the two decay modes are complementary to each other in probing the stop 
pair production, and the LHC with y/s = 14 TeV and 100 fb^ 1 luminosity is capable of discovering 
the stop predicted in natural MSSM up to 450 GeV. If no excess events were observed at the LHC, 
the 95% C.L. exclusion limits of the stop masses can reach around 537 GeV. 

PACS numbers: 14.80.Da,14.80.Ly,12.60.Jv 
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I. INTRODUCTION 



Although the standard model (SM) has been successful in describing the existing exper- 
imental data, it is suffering from the hierarchy problem and new physics based on certain 
symmetry is widely expected to appear at TeV scale to stabilize the electroweak hierarchy 
against radiative corrections. This belief was further strengthened by the recent discovery 
of the Higgs boson at the Large Hadron Collider (LHC) with its mass determined around 
125 GeV [l|. This mass value agrees well with the prediction of low energy supersymmetry 
(SUSY), which is so far the most promising new physics candidate. 

In SUSY, all known bosons and fermions have their supersymmetric partners, and the 
scalar top quarks (called stop £j with i = 1,2), as the top quark partners, can modify the 
property of the SM Higgs boson by exactly canceling out the dangerous quadratic divergence 
of the top quark loop. Obviously, the experimental determination of the stop properties is 
crucial to unravel the nature of supersymmetry in protecting the Higgs mass at the weak 
scale and thus solving the hierarchy problem. In fact, such activities have been carried out 



extensively at the hadronic colliders such as the LHC and the Tevatron 



2|4iJ , but in contrast 



with the strong mass bounds (about 1 TeV) on the gluino and the first generation squarks 
0, a relatively light stop (say about 300 GeV) can not be excluded. Nevertheless, it should 
be mentioned that the recently measured Higgs boson mass around 125 GeV may give some 
indications for the stop sector 0-8]. In the popular MSSM with moderate tan/3 and large 
rriA, the Higgs mass is given by jg] 



mi ~ Ml cos 2 2/3 + In — | + - - * - * 1 , (1) 

where v = 174 GeV, X t = A t — \i cot (3 and M$ is the average stop mass scale defined by 
Ms = y/m^m^. This expression indicates that, for the heavier stop £2 around 1 TeV as 
discussed above, the lighter stop t\ must be heavier than about 200 GeV and \X t \ must 
be larger than 1.5 TeV in order to push the Higgs mass up to 125 GeV 8j. About these 
constraints, one should keep in mind that they are independent of the decay modes of t\, 
but on the other hand, they may be greatly weakened if there exists additional contribution 
to the Higgs mass [8|, |9j . 

On the theoretical side, there are good reasons to consider at least one stop significantly 
lighter than other squarks with a mass around several hundred GeV. Firstly, in some pop- 
ular grand unification models, supersymmetry breaking is usually assumed to transmit to 



the visible sector at a certain high energy scale, and then Yukawa contributions to the 
renormalization group evolution tend to reduce stop masses more than other squark masses. 
Secondly, the chiral mixing for certain flavor squarks is proportional to the mass of the cor- 
responding quark, and is therefore more sizable for stops. Such a mixing will further reduce 
the mass of the lighter stop. Thirdly, in the MSSM the minimization conditions of its Higgs 



potential imply 10 1 



Ml _ {m\ d + £ d ) - (m 2 Hu + S M ) tan 2 g 
2 ~ tan 2 p - 1 



/A (2) 



where rn 2 H and rn 2 Hu represent the weak scale soft SUSY breaking masses of the Higgs 
fields, fi is the higgsino mass parameter, tan (3 = v u / Vj. S u and arise from the radiative 
corrections to the Higgs potential and the dominant contribution to the S u is given by 



31? 
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These two equations indicate that, if the individual terms on the right hand side of Eq. (J2]) 
are comparable in magnitude so that the observed value of Mz is obtained without resorting 
to large cancelations, the natural values of /i and nip should be around 100 GeV and several 
hundred GeV respectively. Numerically, the requirement of E u < M§/2 (or S u < v 2 ) leads 



to Tni 12 upper bounded by about 500 GeV (or 1.5 TeV) [ll|. Moreover, we note that a 



light stop is also phenomenologically needed by the electroweak baryogenesis 



12] and may 



be welcomed by the dark matter physics 13]. In the MSSM, although the gluino contribute 



to m 2 at the one-loop level and to rn 2 Hu at two-loop level, the corrections are proportional 



to m~ and can be greatly enhanced by the large gluino mass 14j. In order to the keep the 
naturalness, we expect m g to be lighter than about 3 TeV for mi < 1.5 TeV. However, since 
the current results of searching for the supersymmetry indicate that a gluino with mass 
about 1 TeV can safely avoid the LHC constraints, we require 1 TeV < nig < 3 TeV in our 
calculation. 

Motivated by the theoretical preference and the results of the LHC search for SUSY, 



recently the natural MSSM scenario attracted broad attention 14j-|l8], which focuses on the 
following parameter space of the MSSM 

• M ^ 100 - 200 GeV and m h 2 < 1 - 1.5 TeV as preferred by Eq.© and Eq.©; 
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• 1 TeV < nig < 3-4 TeV to escape the LHC constraint and at same time to avoid 
spoiling color symmetry; while the electroweak-ino masses may still be at sub- TeV 
scale; 

i 

• rriA ~ \ m H d \ 2 ~ I A* I t an /3 as suggested by the relation m\ ~ 2u 2 + m 2 Hu + m 2 Hd + T lu + T ld 
and Eq.([2]); 

• rriq 12 , a ~ 10 — 50 TeV to provide a decoupling solution to the SUSY flavor and 
CP problems. 

Since the stop are relatively light and sensitive to probing this scenario, there have been 
recently many theoretical studies on the collider signatures of the light directly produced 
stop in the i?-parity conserving and violating MSSm[3]. For example, by using the top 
tagging technique, the sensitivity of stop searches were studied in the hadronic, semi-leptonic 
and di-leptonic channels 



20H23| . In order to suppress the di-leptonic top backgrounds, the 



authors in Ref. 24| explore some new kinematic observables developed from Mt 2 to improve 
the sensitivity of the stop searches. For the small mass splitting between stop and top, it is 
pointed that the rapidity difference and spin correlation of the daughter products from stops 
decay can be helpful to discover the signal{25]. When the stop mass is close to the lightest 
supersymmetric particle mass, the monojet signature from t\t\ + j production is expected 
to be useful in detecting the stop[26j. If the stop mass is degenerate with the sum of the 



masses of its decay products, the searches based on missing transverse energy (E™ tss or 
have significant reach for stop masses above 175 Gevj^]. When the .R-parity is violated, 
the decay modes of the stop will be very different from the ones in i?-parity conserving 
MSSM, such as stop decaying to dilepton and trilepton final states 28[]. We also noted that 
the constrains on the light stop in the natural SUSY have been discussed by using the results 
from sparticles searches at the LHC, and indicated that they were mild and can be safely 



avoided current ly[14j, |17|, 
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29|. 



In this work, we investigate the potential of the LHC in probing the lighter stop t\ 
predicted by the natural MSSM with .R-parity, which is based on some considerations: (i) 
Most of the studies of the stop searches have been carried out under some assumptions 
at the LHC in a model independent way or in simplified models. It will be meaningful 
to explore what might happen in a realistic model like MSSM under the current available 
experimental constraints; (ii) Due to the i?-parity conservation, there will be sizeable missing 



energy appearing in the sparticles productions and decays, which can be easily identified 
in the LHC data; (iii) One interesting phenomenological feature of the natural MSSM with 
-R-parity is that both the lightest neutralino and the lighter chargino are Higgsino-like, and 
consequently t\ always decays dominantly into t\\ and bxf with xt — > XiW*, which can 
greatly simplify the analysis of the ti detection at the LHC. 

For this purpose, we first scan the parameter space of the natural MSSM by considering 
various constraints in Sec. II. Then in Sec. Ill we discuss the observability of t\ through the 
direct stop pair production in the allowed parameter space by performing the Monte Carlo 
simulations for the channel pp — > tit\ — > tXi^X? an d the channel pp — > tit* — > bx\bXi- We 
will present their corresponding sensitivities for 8 TeV LHC and for 14 TeV LHC respectively. 
Finally in Sec. IV, we summarize the conclusions obtained in this work. 
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FIG. 1: Scatter plot of surviving samples in natural MSSM, projected on the plane of mi versus 
\At\/Ms- All the samples here predict the SM-like Higgs boson in the mass range 125 ± 2 GeV. 

II. SCAN OVER THE PARAMETER SPACE 

Motivated by the natural MSSM, we scan the parameter space of the MSSM as follows: 

l<tan/3< 60, 100 GeV < /i < 200 GeV, 100 GeV < M 2 < 1 TeV, 

\A t \ < 3 TeV, 100 GeV < (Mq 3 , M^) < 2 TeV, 90 GeV < M A < 1 TeV. (4) 
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FIG. 2: Same as Fig.l, but showing the relic density of the neutralino dark matter. 
For other unimportant parameters, we fix all the soft breaking parameters in the slepton 



sector and the first two generation sector at 10 TeV, and we assume A t = Aj,, Mq 



M 



and Mi : M 2 = 1:2 (inspired by the grand unification relation). In our scan, we consider 
following constraints: 

• We require the SM-like Higgs mass within the range 125 ± 2 GeV. We use the code 
FeynHiggs2.8.6 |30[ to calculate the mass and the code HiggsBounds-3.8.0 [31] to con- 
sider the experimental constraints on the Higgs sector of the natural MSSM. 



Since the natural MSSM has important implications in B-physics 32], we use the 
code susy_flavor v2.0 [33] to consider the constraints from the processes B — > X s 7 and 



B, 



We consider indirect constraints from the precision electroweak observables such as pi, 
i 



sin 2 l e ff, m w and i4- We use our own code for such calculation (34 ] . 



We require the thermal relic density of the lightest neutralino (as the dark matter 



35j. We use the code MicrOmega v2.4 [36( to 



candidate) is below the WMAP value 
calculate the density. 

After analyzing the surviving samples, we find they have two characters. One is that the 
Higgs mass of 125 ± 2 GeV requires mf > 220 GeV and there is a rather strong correlation 
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FIG. 3: Same as FigfTj but showing the cross section of t\t\ production versus the stop mass, 
between m~ tl and the ratio \A t \/Ms, as shown in FigJTJ Here we further clarify that, if Mq 



The other feature is that for most cases, the values of /z are significantly smaller than M\ 
so that the lightest neutralino is higgsino-like. Fig. 2 indicates that the surviving samples 
lie within two isolated regions. We checked that the lightest neutralino is bino-like in the 
left region and higgsino-like in the right region. Here the bino(higgsino)-like means it is 
still a mixed state but the dominant component is bino(higgsino). For the light neutralino 
dark matter (bino-like), the main annihilation channel is through exchanging Z boson. The 
annihilation cross section is roughly proportional to l/(4m^ — m 2 z ) 2 . When the neutralino 
mass is about 50GeV— 60GeV, the annihilation cross section may be very large, so the 
relic density will be less than 0.1. When the neutralino becomes heavier (60GeV—90GeV, 
neutralino is still bino-like), the annihilation cross section will drop. The relic density 
becomes large and even exceeds the WMAP value, and these samples are excluded. This 
is the reason for the gap between 60GeV— 90GeV. When the neutralino goes on becoming 
heavy(>90 GeV), the dominant component of the neutralino will be higgsino. The coupling 
between neutralino and Higgs gets important and annihilation rate goes up, then the relic 
density drops. 

About the natural MSSM, we have two comments. One is that in this scenario the di- 
photon signal of the SM-like Higgs boson can hardly be enhanced to satisfy the requirement 
of the LHC data. This is because in the framework of the MSSM, there are only two cases 
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which can enhance the di-photon rate, i.e. the small a e ff scenari o |37l . |38| and the light f 



scenario [6|, I39f . and in each case a large \x is needed. In the Ref. |40j], the authors pointed 
that the light stop with large couplings to Higgs boson in the SM+ stop model can improve 
the SM fitting to the LHC and Tevatron data by enhancing F(h — > 77) and suppressing 
r(/i — > gg). However, we should note that it does not mean that the di-photon production 
rate can reach the measurement of the LHC in a concrete MSSM model since the reduction 
of T(h —7- gg) is usually much stronger than the enhancement of T(h — > 77) for large values 



of X t (X t = A t — /itan/3) 41] . The other is that recently the ATLAS collaboration searched 
for the gluino- mediated stop pair production followed by the decay ti — > bxt — > b£i/x®, 
which set a lower bound > 450 GeV [3]. This conclusion is not applicable to the our 
calculations since we take the gluino mass to be larger than ITeV in the allowed parameters 
space of the natural SUSY. 

III. OBSERVABILITY OF STOP PAIR PRODUCTION AT THE LHC 

In this section we discuss the LHC potential of discovering the stop through the direct 
stop pair production in the natural MSSM at -y/i = 8, 14 TeV. In Figj3]we show the pp — > tit* 
production rate at the next leading order for the surviving samples. In getting this figure we 



used the package Prospino2.1 |42|] and the parton distribution function CTEQ6.6m |43] with 
the renormalization scale fiR and factorization scale fiF setting to mf . This figure indicates 
that the maximal values of the cross section reach 5.5 pb and 25.7 pb for the LHC with 
yfs = 8 TeV and y/s = 14 TeV respectively, and with the increase of the stop mass, the 
production rates drop rapidly. 

In FigJH we present various decay branching ratios of ti which are obtained by using 



the package SDECAY 44j. This figure indicates that for m~ tl < 320 GeV where the decay 
channel t\ — > t\\ does not open up, t\ decays into bxt whh a ratio of 100%, and as the stop 
becomes heavier, the branching ratios for ti — > bxt an d t\ — > txl may still be around 50%. 
In contrast, the branching ratios for t\ decays into t%g 4 and bxt are usually less than 20%. 

In the following we perform detailed Monte Carlo simulations to investigate the observ- 
ability of the direct stop pair production at the LHC. We concentrate on the semi-leptonic 
analysis with the b-tagging efficiency 40%, where the signal is consisted of four jets(at least 
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FIG. 4: Same as Figj3j but showing the decay branching ratios of the stop, 
one b-jet), one lepton (e or /i), and missing transverse energy. We first consider the process 
pp t x t\ (txlml) -> (b£ + vXi)(bjjXi) or (&#X?)(&ri/x°). (5) 

From the ATLAS search for the signal ti + E™ lss j^j , we can see that the dominant SM back- 
ground after the E™ lss and cuts is ti di-leptonic channel with one lost lepton and two 
additional jets from initial state radiation to fake the hadronic W. Another backgrounds in- 
clude tt semi-leptonic channel, tt di-leptonic channel with one r from top decay misidentified 
as a jet, PF+jets and tiZ. Here we emphasize that the tiZ background becomes important 
for a heavy stop and should be considered in estimating the significance. In our calculation, 
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to 



we normalize the signal and the ti back grou nd to their NLO values J42l. |45l| . and simulate the 
signal and backgrounds by MadGraph5 j46| interfaced with PYTHIA 47] and Delphes 
carry out the parton shower and fast detector simulation. We use the anti-fc* algorithm 
with the distance parameter R = 0.4 to cluster jets and the MLM scheme 50[ to match our 
matrix element with parton shower. We checked that the shapes of the matched VF+1,2,3 
partons are very similar, and for simplicity, we take W+2 jets samples in our calculations. 
In our calculations, since we employ the the variable M^f 2 defined in Ref . [24| , we checked our 
results with theirs for the same parameters at yj~s = 7 TeV and found they were consistent 
with each other. 



In Figj5l we show the distributions of $ T , the transverse mass Mt defined 
Mjf 2 for the backgrounds and our benchmark point m+ 
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FIG. 5: The distributions of £+4-jets+j£'y with respect to $Ti Mt and Mjf 2 for the signal pp — > 
t\t* — > tXitXx and backgrounds. In the Mjf 2 distribution, we impose the cuts $ T > 150 GeV and 
Mt > 150 GeV on the events and only display the di-leptonic tt, tiZ backgrounds and signal, 
where the events with wrong or no solution for are included in the last bin. 



TABLE I: The significance of stop pair production pp — > tit* — > txitXi f° r 100 fb 1 luminosity 
after imposing various cuts. Here we take = 429 GeV and x? = HO GeV for illustration. 
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with = 8 TeV (similar results are found for yfs = 14 TeV). This figure indicates that 
most events of W+jj and semi-leptonic tt backgrounds are characterized by $t — 100 GeV 
and Mt < 100 GeV, and most events of the di-leptonic tt backgrounds are characterized 
by Mj^ < 170GeV, while a significant fraction of the signal may have larger $ T) Mt and 
M^. Figj5]also indicates that the distributions of the ttZ background are quite similar to 
the signal and are difficult to be suppressed. Fortunately, the production rate of tiZ is much 
smaller than the one of tit\. 

In Table I, we present the significance S/ \f~B of our benchmark point for 100 fb~ x lumi- 
nosity with ^/s = 8 TeV and 14 TeV respectively by sequentially imposing the cuts on $ T , 
Mt and M^. It can be seen that, for the given reference point, the cut Mt > 150 GeV can 
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greatly enhance the significance and M^ 2 > 173 GeV further improves the significance by 
about 14% for y/s = 8 TeV and 9% for = 14 TeV to reach 3.11 and 8.60 respectively. 

Therefore, for our simulations in the allowed parameters space, we take the following 
events selection criteria: 



• One isolated electron or muon that passes the following requirements; 

- Electrons E T > 25GeV and \r]\ < 2.47 without 1.37 < \r]\ < 1.52; 

- Muon: E T > 20GeV and \r]\ < 2.5; 

- Events are rejected if they contain a second lepton candidate with P? > 15GeV; 

• Four or more reconstructed jets with Pt > 25GeV and \t)\ < 2.5. 

• $ T >150GeV, M T >150GeV, >173GeV. 

where the basic cuts about pt and rj on leptons and jets are from the ATLAS report In 
order to improve the signal sensitivity, we increase the values of ATLAS cut $ T from 100 
GeV to 150 GeV to further suppress the semi-leptonic ti background and use the new cut 
Mjf 2 >173 GeV to reduce the di-leptonic ti background in our calculations. 
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FIG. 6: The significance of stop pair production pp — > t\t* — > tx\ix\ for the surviving samples in 
the natural MSSM. 



In Fig. 6 we show the significance of the surviving samples with < 800 GeV. This 
figure indicates that the largest significance can be reached at m t - ~ 430 GeV where the 
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FIG. 7: Same as Figj5j except that in plotting Mjf 2 distribution we impose the cuts $ T > 100 
GeV and M T > 100 GeV on the events. 

significance is about 1.5 for y/s = 8 TeV with 20 fb~ l luminosity and 8.5 for y/s = 14 TeV 
with 100 fb -1 luminosity, and with the increase of the stop mass, the significance drop by 
one half for ~ 500 GeV mainly due to the reduction of the production rate. Our results 



are not as optimistic as those in [20l . |21 



24j because we have taken into account the branch 
ratio of t\ — > t\\. Figj6]also indicates that there are two branches for the significance in the 
mass region 320 GeV< < 600 GeV. We checked that the upper branch corresponds to 
high branching ratio of t\ — >■ txl, which varies from 42.1% to 57.7% and results in a large 
signal rate, while the lower branch corresponds to a small ratio due to the competition of 
the decay mode t\ — > txl with t\ — > bxf and ti —> tx®- 

Above analysis implies that, in order to fully explore the parameter space of the natural 
MSSM in stop detection, the decay mode t\ — > bxl should also be considered. So we next 
consider the process 

pp -> U\ -> (bxt)(bxi) -)• (b£ + PX i)(bjjfi) or (bjjxi 



(6) 



Same as in Figj5j we show the distributions of the three variables in Fig J7] for the benchmark 
point = 273.6 GeV, m-+ = 163.5 GeV and m^o = 156.3 GeV. Compared with the 
distribution in FigfSl one can see that more signal events have lower values of p T and 
low M T , due to the relatively light t\. FigfT] also indicates the M™ 2 variable is helpless in 
suppressing the di-leptonic ti events and any cut on Mj^ may hurt the signal greatly. 
In FigJHl we show the significance of the surviving samples for the process in Eq.([6]). In 
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order to keep more signal events, here we relax the cuts of and Mt used for the process 
in Eq.flH]) as follows: 

$ T > 100 GeV, M T > 100 GeV. (7) 

From this figure, one can learn that, due to the large stop pair production rate, the signifi- 
cance for mi x = 250 GeV may reach 7 for yfs = 8 TeV with 20 fb^ 1 luminosity and 64 for 
y/s = 14 TeV with 100 fb~ l luminosity, but for m~ tl = 400 GeV, the maximum value drops 
to 1.5 and 10 respectively. 
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FIG. 8: The significance of stop pair production pp — > t±t* — > (bxf)(bxi ) for the surviving samples 
in the natural MSSM. 



Finally, we summary the significance of the direct stop pair production with the above 
two decay modes of t\ for y/s = 14GeV and 100 fb~ x luminosity. The results are displayed 
in Figj9] where only the maximal significance under each cut is shown. This figure indicates 
that, for <400GeV, detecting the stop pair production through the chargino decay is 
more effective, while for 400GeV < m~ tl < 450GeV the neutralino decay is more effective. 
This figure also indicates that the LHC can discover t% predicated in nature MSSM up to 
450GeV. If no excess events were observed at the LHC, the 95% C.L. exclusion limits of the 
stop masses can go up to around 537 GeV no matter what decay modes of the stop in the 
natural MSSM. 
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FIG. 9: The maximal significance in Eq.© or Eq.© as a function of the stop mass. The bullets 
(red) are obtained in Eq.© while the triangles (green) are obtained in Eq.([5]). 

IV. CONCLUSION 

In this work we studied the direct stop pair production at the LHC in the natural MSSM. 
We scanned over the corresponding parameter space by considering various experimental 
constraints and then in the allowed parameter space we examined the observability of the 
direct stop pair production at the LHC through the semi-leptonic analysis. We focused on 
the following two channels 

pp tJl (*X?)(tx?) ->• (bl + vXi)(bjjX°i) or (&TPx°), 
pp tJl (bxtKbxi) (b£ + ^Xi)(bjjX i) or (6jjx?)(&r *x°), 

and performed detailed Monte Carlo simulations about the signals and backgrounds. We 
found that for < 400GeV the second channel is better while for 400GeV < mi < 
450GeV the first channel is better. We also found that the LHC with y/s = 14TeV and 
luminosity is capable of discovering t\ predicated in nature MSSM up to 450GeV. 
If no excess events were observed at the LHC, the 95% C.L. exclusion limits of the stop 
masses can reach around 537 GeV in the natural MSSM. 
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FIG. 10: The ATLAS observed and our expected exclusion limit in the plane of m^o versus m% 1 , 
assuming BR{h txi) =100%. 

Note Added 

Very recently, the ATLAS collaboration reported the result of the direct searching for the 
stop pair production base on 4.7/6 -1 of data[5l|. We validated our simulation by reproducing 
the ATLAS exclusion limit according to the assumptions and cuts in the report as follows: 

• One isolated electron or muon passing 'tight' selection criteria; 

- Electrons E T > 25GeV and \rj\ < 2.47; 

- Muon: E T > 20GeV and \r]\ < 2.4. 

• Four or more jets with |7/| < 2.5 and Pt >80, 60, 40 and 25 GeV, and at least one jet 
to be identified as a b-jet; 

• A(j) min > 0.8,where A0 min is the minimum azimuthal separation between the two 
highest Pt jets and the missing transverse momentum direction; 
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• The jet-jet pair having invariant mass > 60 GeV and the smallest AR is selected to 
form the hadronically decaying W boson. The mass rrijjj is reconstructed including 
a third jet closest in AR to the hadronic W boson momentum vector and 130 GeV 
< rrijjj < 205 GeV is required; 

• E™ iss >150 GeV, E™ iss /^H^ >7 GeV 1 / 2 and m T > 120 GeV; 

• Events are rejected if they contain additional leptons passing looser selection criteria. 
Here we treat the looser selection criteria as P T > 15 GeV. 

• The branch ratio of ti — > tx® is assumed to be 100%. 

In Fig. 10, we display the expected exclusion limit from our simulation. Considering the 
differences between the fast simulation and full detector simulation, we can see that our 
result is consistent with the ATLAS exclusion limit within the reasonable error range. We 
also expect our result can be improved by the simultaneous fits method used by ATLAS 
for five signal regions and three control regions, however, which is beyond the scope of our 
simulation. It should be noted that the above stop masses limits can be avoided in our 
study, since the stop decays with a mixture of the branching ratios. 
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